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When a quantum many-particle system exists on a randomly diluted lattice, its intrinsic 
ly^ ' thermal and quantum fluctuations coexist with geometric fluctuations due to percola- 

tion. In this paper, we explore how the interplay of these fluctuations influences the 
■ phase transition at the percolation threshold. While it is well known that thermal fluctu- 

ations generically destroy long-range order on the critical percolation cluster, the efli'ects 
of quantum fluctuations are more subtle. In diluted quantum magnets with and without 
, dissipation, this leads to novel universality classes for the zero-temperature percolation 

(-H 1 quantum phase transition. Observables involving dynamical correlations display nonclas- 

Q ' sical scaling behavior that can nonetheless be determined exactly in two dimensions. 
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1. Introduction 



> 
OO 

' In disordered quantum many-particle systems, random fluctuations due to impuri- 

ties and defects coexist with quantum fluctuations and thermal fluctuations. Close 
^ to phase transitions, the interplay between these different types of fluctuations 

■ can cause many unconventional phenomena such as quantum Grifhths effects, ^'^ 

I infinite-randomness critical points'^'* and the destruction of the phase transition by 

' smearing^ (for a recent review see, e.g., Ref. 6). A particular interesting case of this 

scenario are randomly diluted magnets. Site or bond dilution defines a percolation 
problem'' for the lattice which can undergo a geometric phase transition between a 
^ • disconnected and a percolating phase. 

5^ , Here, we discuss how the interplay between the quantum fluctuations of the 

spins and the geometric fluctuations of the lattice changes the phase diagram and 
the nature of the magnetic percolation transition. Our paper is organized as follows: 
In Sec. 2, we collect the basic results of classical percolation theory^ to the extent 
necessary for the following sections. Section 3 is devoted to the behavior of classical 
magnets on percolating lattices. These are older results summarized here mainly 
for comparison with the quantum case. Section 4 is the main part of the paper. 
Here, we discuss three different examples of percolation quantum phase transitions 
(QPTs) of diluted quantum magnets. We conclude in Sec. 5. 
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Fig. 1. Site-diluted square lattice at impurity concentrations below (p = 0.3), at {p = Pc ^ 
0.4073), and above (p = 0.6) the percolation threshold. 

2. Geometric Percolation 

Classical percolation theory'' deals with the geonicitrie properties of randomly di- 
luted lattices. The central question is whether the diluted lattice contains a large 
cluster that spans the entire sample or whether it is decomposed into small dis- 
connected pieces. For definiteness, consider a d-dimensional hypercubic lattice with 
bonds between the nearest neighbor sites in which a fraction p of all sites (site per- 
colation) or bonds (bond percolation) is removed at random. Depending on p, the 
lattice can be in one of two "phases" , separated by a sharp percolation threshold at 
p = Pc- li p < Pc (the percolating phase), there is a single large cluster that spans 
the entire sample (as well as some smaller clusters). In the thermodynamic limit, 
this cluster, the so-called infinite cluster, becomes infinitely large and contains a 
nonzero fraction P^o of all sites. In contrast, for p > Pc the lattice is decomposed 
into small disconnected finite-size clusters only (sec Fig. 1). Right at p = Pc, there 
are clusters on all length scales, and their structure is fractal. 

The behavior of the diluted lattice close to Pc is very similar to the critical 
behavior near a continuous (2nd order) phase transition with the geometric fluctua- 
tions due to dilution playing the role of the usual thermal or quantum fluctuations. 
This implies that observables are governed by power-law scaling relations. A central 
quantity is the cluster size distribution ns{p). It measures the number of connected 
clusters with s sites (per lattice site). Close to Pc, it takes the scaling form 

n«(p) = s-^=/[(p-Pe)s"1 , (1) 

where CTc and Tc are critical exponents. In two dimensions, they are known ex- 
actly, CTc = 36/91 and Tc = 187/91; and in three dimensions they are well 
known numerically, ctc « 0.45 and Tc w 2.18. f{x) is a scaling function which 
behaves as f{x) ~ exp{—Bix^/'^^) for a; > 0, f{x) = const for x = 0, and 
f{x) ~ exp[— (^22;^/'^'=)^"^/'*] for a; < 0. The critical behavior of all other geo- 
metric properties can be expressed in terms of the exponents ctc and Tc. In particu- 
lar, in the percolating phase, the fraction of sites in the infinite cluster behaves as 
^cx> ~ {Pc — pY" with the exponent /3c given by /3c = (tc — 2)/ac- When approach- 
ing the percolation threshold, the typical linear size of the finite-size clusters, the 
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Fig. 2. Schematic phase diagrams for classical and quantum magnets on diluted lattices. In the 

classical case, magnetic order is destroyed by increasing the temperature T, in the quantum case 
by increasing the quantum fluctuations, e.g., the transverse field hx in a quantum Ising magnet. 

connectedness length, diverges as ^ \p^Pc\~^'' with Vc = {tc — l)/{dac)- Finally, 
the fractal dimension D / of the infinite cluster at the percolation threshold can be 
expressed as Df = d/{Tc — 1). 

3. Diluted classical magnets 

In this section we briefly summarize the behavior of a classical Ising or Heisenberg 
magnet on a randomly diluted lattice. Consider the Hamiltonian 



where Si is a classical Ising or Heisenberg spin at site i, and J > is the exchange 
interaction between nearest neighbors. The dilution is implemented via quenched 
random variables taking the values and 1 with probabilities p and 1 — p, re- 
spectively. We first discuss the temperature-dilution phase diagram. In the absence 
of dilution, the model orders ferromagnetically below a critical temperature Tc(0) 
(provided d > 2 for Ising and c? > 3 for Heisenberg spins). Upon dilution, mag- 
netic order is weakened, and decreases. An important question is whether the 
magnetic phase is completely destroyed before the dilution reaches Pc, right at Pc, 
or whether long-range order survives even on the critical percolation cluster at Pc, 
corresponding to phase diagrams (a), (b), and (c) in Fig. 2, respectively (long-range 
order cannot survive for p > Pc because the system consists of small disconnected 
clusters). Phase diagram (a) can be excluded because the infinite percolation clus- 
ter is a massive d-dimensional object for any p < pc- Since the critical percolation 
cluster at p = Pc has a fractal dimension £)/>!, one might be tempted to conclude 
that it supports magnetic long-range order, at least in the Ising case (implying a 
phase diagram of type (c)). However, this is incorrect: At the percolation thresh- 
old, thermal fluctuations immediately destroy the magnetic order. ^"^'^ This can be 
understood by considering the "red sites" of the critical percolation cluster, i.e., 
sites that divide the cluster into two otherwise disconnected pieces (see Fig. 3). The 
orientation of the spins on these two pieces can be flipped with respect to each other 
with a finite energy cost of 2J. For any T ^ 0, both the parallel and antiparallel 
configuration at each of the red sites contribute to the statistical sum, destroying 
magnetic long-range order. The phase diagram is thus of type (b), and a percolation 
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Fig. 3. Red sites vs. red lines (shown in black) in a critical percolation cluster. 

transition only occurs at exactly zero temperature. Since there are no thermal fluc- 
tuations, the critical behavior of this transition is identical to geometric percolation. 
Any T ^ Q destroys the percolation critical behavior, instead the transition is in 
the universality class of the corresponding generic disordered classical magnet. 

4. Diluted quantum magnets 

We now turn to the main topic, QPTs on percolating lattices. Generally, QPTs occur 
at zero temperature as functions of pressure, magnetic field or other nonthermal 
control parameters (for reviews, see, e.g., Refs. 11-14). One important aspect of 
these transitions is the so-called quantum-to-classical mapping. It arises because in 
quantum statistical mechanics the partition function does not factorize in potential 
and kinetic parts. Instead, it has to be formulated in terms of space and time- 
dependent variables. As a result, (imaginary) time acts like an additional coordinate, 
and a QPT in d dimensions can be related to a classical transition in a higher 
dimension,'^ a fact we will be using repeatedly below. 

4.1. Transverse-field Ising model 

The first example is a randomly diluted Ising model in a transverse magnetic field, 
given by the Hamiltonian 

Hi = -jJ2 ^^^j -h,j2^^^^ ^^Y. ^^^^ ■ (3) 

and Sf are the z and x components of the the quantum spin-1/2 at site i; hx 
is the transverse field that controls the quantum fluctuations, and H is the field 
conjugate to the order parameter. The clean model {p = 0) is a paradigm for the 



''This mapping is restricted to the thermodynamics only. Moreover some quantum transitions lead 
to extra complications such as Berry phases that do not have a classical counterpart. 
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Fig. 4. Schematic T — hx —p phase diagram of a diluted quantum magnet. There is a multi-critical 
point at (T = 0,h*,pc)- The QPT across the dashed line is the topic of this paper. 

study of QPTs: For <^ J, the ground state is ferroniagnetically ordered in z- 
direction while for hx ^ J the quantum fluctuations due to the transverse field 
destroy the long-range order. The two phases are separated by a QPT at hx ~ J-^^ 

As in the classical case, we ask how the dilution influences the phase diagram 
(in the hx-p plane), i.e., is the phase diagram of type (a), (b), or (c) in Fig. 2? 
As before, (a) can be excluded because the infinite percolation cluster is a mas- 
sive rf-dimensional object for p < pc- To decide between (b) and (c), we adopt 
the "red-site" argument to the quantum case. Following the quantum-to-classical 
mapping, wc have to consider an effective system in d space dimensions and one 
extra imaginary time dimension which becomes infinite for temperature T — > 0. 
Crucially, the defect positions are time-independent. Instead of red sites we thus 
have "red lines" separating different pieces of the critical percolation cluster (right 
panel of Fig. 3). Configurations with different spin orientations on two such pieces 
now come with a infinite effective energy penalty and are suppressed. This suggests 
that magnetic long-range order can survive on the critical percolation cluster if the 
quantum fluctuations are not too strong, implying a phase diagram of type (c). This 
has been confirmed by simulation results not only for quantum Ising models but 
also Heisenberg magnets and quantum rotors. 

Combining the effects of thermal and quantum fluctuations, wc obtain the T — 
hx— P phase digram shown in Fig. 4. A diluted quantum magnet can thus undergo 
two nontrivial QPTs, separated by a multicritical point at (T = 0, h*,Pc). We are 
interested in the percolation transition at Pc and hx < h*, i.e., the transition across 
the dashed line in Fig. 4. It was first investigated in detail by Senthil and Sachdev:^^ 
Consider a single percolation cluster of s sites. For small transverse field hx < h^, 
all spins on the cluster are parallel. The cluster thus acts as a two-level system with 
an energy gap (inverse susceptibility) that depends exponentially on the cluster 
size As ~ x7^ ~ hxexp{—Bs) with B \n.{J/hx)- (AH other excitations have at 
least energy J.) Since the cluster size s and its linear extension L are related via 
the fractal dimension s ~ L'-'f , we obtain an unusual exponential relation between 
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length and (inverse) time scales which is sometimes called activated scaling, 

HK/^s) ~ L''^ . (4) 

The critical behavior of the total system can now be found by summing over all 
percolation clusters via the cluster size distribution (1). Let us first consider static 
quantities like magnetization or magnetic spatial correlation length. For hx < h*, 
magnetic long-range order survives on the infinite percolation cluster, while all 
finite-size clusters do not contribute. Thus, the total magnetization is proportional 
to the number of sites in the infinite cluster, m ~ Poo ~ {Pc — p}^' for p < pc- 
The magnetic order parameter exponent /? is thus identical to that of geometric 
percolation. A similar argument can be made for the magnetic correlation length 
^: For hx < h*^, all spins on a cluster are correlated, but the correlations cannot 
extend beyond the cluster size, thus ^ ~ ^ |p — Pc|~''°5 and the magnetic correla- 
tion length exponent is identical to the geometric one, too. In contrast, quantities 
involving quantum dynamics behave nonclassically. For instance, the dependence of 
the magnetization on the ordering field H takes the scaling form 

m(p -pc,H) = b-'^^'"^ m. {[p- pc)h^^'''' , ln(i? ) , (5) 

with h being an arbitrary scale factor. At the percolation threshold p = Pc, this 
gives the unconventional relation m ~ [ln(i?)]^~'^'=. For p ^ Pc, the transition is 
accompanied by strong power-law quantum GriSiths efiFects.^'^^ 

4.2. Bilayer quantum Heisenberg magnet 

This subsection is devoted to diluted Heisenberg magnets. Specifically, we consider 
a dimer-diluted bilayer quantum Heisenberg antiferromagnet with the Hamiltonian 

Hh = J\\ ^ CiejSi^a ■ ^j,a + J± ^ eiSj^i • Sj,2 , (6) 

{i,3) i 
a=l,2 

where Sj^a is the spin operator at site j in layer o = 1 or 2. The clean system 
{p = 0) undergoes a QPT between a paramagnetic and an antiferromagnetic phase 
as a function of the ratio J±/J\\ between the inter-layer coupling and the in-plane 
interaction. For J± ^ Jy, the corresponding spins in the two layers form a singlet 
which is magnetically inert. Thus, there is no long-range order. In contrast, for 
J|| ^ J±_ , each layer orders antifcrromagnetically, and the weak inter-layer coupling 
establishes antiferromagnetic order between the layers. The phase diagram of the 
dimer-diluted system has been determined by Sandvik^^ and Vajk and Greven;^^ it 
is shown in Fig. 5. In agreement with the general arguments given in the last sub- 
section, long-range order can survive at p = pc giving rise to a nontrivial percolation 
QPT (across the short vertical line in Fig. 5). To study this transition, we first map 
the low-energy physics of (6) onto a quantum rotor model with the action^^ 

■A= I (ir^ J||eiejSi(r) • Sj(T) — ^^eia;^Si(u;„)Si(-u;„) . (7) 
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Fig. 5. Phase diagram of the dimer-dilutcd bilayer Heisenberg antiferromagnet (after Refs. 16,19). 
Inset: Sketch of the system. 

Here, each rotor variable Si(T) (a unit vector at site i and imaginary time r), de- 
scribes a dimer of corresponding spins in the two layers. a;„ is a Matsubara frequency, 
and the parameter g is related to the ratio Ja_/J\\ of the quantum Hamiltonian (6). 

Our approach^" is the same as in Sec. 4.1, we first consider a single percolation 
cluster of size s and then sum over all clusters by means of the cluster size distribu- 
tion (1). For small g, all rotors on the cluster are correlated but collectively fluctuate 
in time. Thus, each cluster acts as a (O-l-l)-dimensional rotor model with magnetic 
moment s. Its low-energy properties can be easily found by a renormalization group 
calculation or dimensional analysis, leading to a scaling form of the free energy 

{g, H, T) - (.g/.)<D (i/sVs, Ts/g) (8) 

as a function of g, T, and magnetic field H. Here, $ is a universal scaling function. 
Eq. (8) implies that the thermodynamics of a quantum spin cluster is more singular 
in its size s than that of a classical cluster. In particular, classically, the magnetic 
susceptibility increases like x's ^ ■^^ while in our quantum model at T = 0, it 
increases more strongly, Xs ^ . The dynamical critical exponent can be obtained 
by relating the gap A to the linear size L of the cluster via Xs ~ s^/A giving 
A ~ s^^ ^ L^^f . Thus, the dynamical exponent at the percolation QPT is z = Dj. 
The total free energy is obtained by summing (8) over all clusters. This gives rise 
to the general scaling scenario 



2-a = + , 


(9) 


p={d-Df)iy, 


(10) 


7 = {2Df-d + z)v , 


(11) 


S={Df + z)/{d-Df) , 


(12) 


2-1] ^ 2Df -d + z . 


(13) 



All critical exponents are completely determined by two geometric percolation expo- 
nents (say Df and v = Vc) together with the dynamical exponent z which contains 
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Table 1 . Critical exponents of the geometric and quantum 
percolation transition in two and three dimensions.^" 







2d 




3d 


exponent 


classical 


quantum 


classical 


quantum 


a 


-2/3 


-115/36 


-0.62 


-2.83 


P 


5/36 


5/36 


0.417 


0.417 


7 


43/18 


59/12 


1.79 


4.02 


5 


91/5 


182/5 


5.38 


10.76 


V 


4/3 


4/3 


0.875 


0.875 


V 


5/24 


-27/16 


-0.06 


-2.59 


z 




91/48 




2.53 



the information on the quantum fluctuations. Thus, a, 7, (5, and r\ are modified 
by the quantum dynamics while /3 and v are unchanged. The resulting exponent 
values are shown in table 1. We have recently confirmed the 2d results by per- 
forming large-scale Monte-Carlo simulations of the action (7).^^'^^ Let us note that 
the behavior of sitc-dilutcd (rather than dimcr-dilutcd) Hcisenbcrg antifcrromag- 
nets is more complicated because the effective action contains Berry phases. Recent 
computer simulations^^ suggest that z « l-SD/ to 2Z)/ in this case. 



4.3. Percolation and dissipation 

In many real systems, magnetic degrees of freedom are coupled to a dissipative 
environment of "heat bath modes" (e.g., electronic degrees of freedom in a metal 
or nuclear spins in molecular magnet). In this subsection, we study the influence of 
Ohmic dissipation on a percolation QPT.^^ To this end we add baths of harmonic 
oscillators to the diluted transverse-field Ising model of Sec. 4.1, 



(14) 



where ai,„ and a\ ^ are the annihilation and creation operators of the n-th oscillator 
coupled to spin i; Vi^n is its frequency, and Ai,„ is the coupling constant. All baths 
have the same spectral function £{ijj) = ttX^^ A?„5(a; — Vi,n)/'^i,n = 2nacje~'^/'^'' 
with a the dimensionless dissipation strength and cJc the cutoff energy. 

Following our general approach we first consider a single percolation cluster of 
size s. Without dissipation, its low-energy properties are described by a quantum- 
mechanical two-level system (see Sec. 4.1). In the presence of the heat baths, the 
cluster therefore behaves as a dissipative two-level system with effective dissipation 
strength sa. The physics of this problem is very rich, it has been reviewed, e.g., 
in Ref. 24. For our purposes, the most important aspect is that with increasing 
dissipation strength, the (Ohmic) dissipative two-level system undergoes a phase 
transition from a fluctuating phase at sa < 1 to a localized (frozen) phase at 
sa > 1. As a result, for any given microscopic dissipation strength a, the total 
diluted lattice consists of a mixture of large frozen clusters that act as classical 
superspins and smaller clusters that behave quantum mechanically down to the 
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Fig. 6. Schematic ground state phase diagrams of the diluted transverse-field Ising magnet with- 
out (a) and with (b) dissipation. CSPM is the classical superparamagnetic phase (after Ref. 23). 

lowest temperatures. The resulting phase diagram^^ of the of the dissipative diluted 
transversc-ficld Ising model is shown in Fig. 6. 

The behavior of observables close to the percolation transition can be obtained 
by summing the results for the dissipative two-level system over the cluster size dis- 
tribution (1). The total magnetization has 3 parts: The infinite percolation cluster, 
if any, contributes moo ~ -foo ~ {Pc —pY"- The frozen finite size clusters indi- 
vidually have nonzero magnetization, but they do not align in the absence of an 
ordering field. Finally, the small fluctuating clusters have vanishing magnetization. 
The interplay between these 3 contributions and an ordering field leads to exotic 
hysteresis phenomena. ■^'^ The low-temperature susceptibility is dominated by the 
frozen clusters and behaves classically, X \P ~ PcC" /T. In contrast, the specific 
heat is determined by quantum fluctuations giving C ~ l/\T?{hx/T). 

5. Conclusions 

In summary, we have discussed the interplay of geometric, thermal and quantum 
fluctuations at the percolation threshold. While thermal fluctuations immediately 
destroy magnetic long-range order on the critical percolation cluster, the effects of 
quantum fluctuations are more subtle. Generically, magnetic long-range order on 
the critical percolation cluster can survive a finite amount of quantum fluctuations. 
This gives rise to a nontrivial percolation QPT and a multicritical point separating 
it from the generic "disordered" transition a± p < pc- 

We have discussed three examples of such percolation QPTs in quantum Ising 
and Heisenberg magnets with and without dissipation. In all cases, the critical 
behavior is different from classical (geometric) percolation, but it can be expressed 
in terms of the geometric percolation critical exponents. Percolation transitions are 
thus among the very few examples of QPTs with exactly known exponent values 
in two dimensions. This is caused by the fact that at our percolation transitions, 
the criticality is due to the geometric criticality of the underlying diluted lattice. 
However, the quantum fluctuations "go along for the ride" and modify the behavior 
of all quantities involving dynamic correlations. 

In the quantum Heisenberg case (Sec. 4.2), this leads to new critical expo- 
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nents while the overall power-law scaling scenario is still valid. In contrast, for 
the transverse- field Ising model (Sec. 4.1), the dynamical scaling is of activated 
(exponential) rather than the usual power-law type. Finally, in the presence of dis- 
sipation (Sec. 4.3), the singularities in the dynamics become even stronger such that 
individual finite-size clusters can undergo the phase transition independently from 
the bulk. This leads to a novel superparamagnetic classical cluster phase. Note that 
these different cases agree with a general classification of phase transitions in the 
presence of disorder^ based on the effective dimensionality of the defects. 
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